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An heuristic semiclassical procedure that incorporates quantum gravity induced correc- 
tions in the description of photons and spin 1/2 fermions is reviewed. Such modifications 
are calculated in the framework of loop quantum gravity and they arise from the granular 
structure of space at short distances. The resulting effective theories are described by 
power counting non-renormalizable actions which exhibit Lorentz violations at Planck 
length scale. The modified Maxwell and Dirac equations lead to corrections of the en- 
ergy momentum relations for the corresponding particle at such scale. An action for the 
relativistic point particle exhibiting such modified dispersion relations is constructed and 
the first steps towards the study of a consistent coupling between these effective theories 
are presented. 

1. Introduction 

Recently there has been a revival in the interest of studying bothiihservational and 
theoretical manifestations of quantum gravity induced effects. SoBa This means to 
consider phenomena associated with the Planck scale, which are highly suppressed 
in standard scenarios. Many of the envisioned effective theories for particles which 
incorporate Planck scale modifications also induce minute violations of Lorentz 
covariance at such scale. In this way, these studies naturally overlap with the 
systematic approach developed by Colladay and Kostelecky u which provides the 
most general power counting renormalizable extension of the standard model that 
incorporates both Lorentz and CPT violations. This framework has been used to 
set experimental bounds upon the interactions that produce such violations an; 
the observations performed so far cover a wide range of experimental settings 
As pointed out some time ago, the propagation of high energy particles through 
cosmological distances may provide a realistic possibility to observe one of the effects 
associated with Planck scale corrections: the modification of the corresponding 
energy-momentum relations.13 In the case of photons the authors of Ref.(l) proposed 
to consider the modified dispersion relations 

c 2 p 2 = E 2 {1 + ZE/E QG + 0(E I Eqg) 2 ), (1) 
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with Eqc being a scale expected to be close to the Planck mass Mp = 1/ip — 
10 19 GeV . The relation (Q) implies an energy dependent photon velocity leading to 
a time retardation between two photons simultaneously emitted which are detected 
with a difference in energies AE after traveling a distance L. The uncorrected 
expression for such retardation is 

At « £ (AE / Eqg) [Lj c). (2) 

In the process of detecting photons from the active galaxy Markarian 421 (L = 
3.5 l.y.), Biller et. alE have identified events with AE — ITeV arriving to earth 
within the time resolution of the measurement: At = 280s. In this way the lower 
bound Eqg/£, = 4 x 10 16 GeV is established. Time resolutions £>f milliseconds 
have been achieved in recent gamma ray burst (GRB) observations^ and they will 
substantially improve up to 10~ 7 s with the Gamma Ray Large Area Telescope 
(GLAST) to be operating in the International Space station by 2006. Nevertheless, 
the cosmic radiation background (CRB) prevents space to be transparent to the 
propagation of very high energy photons. For this reason, the detection of ultrahigh 
energy neutrinos (UHEN) could provide an arena to observe such effects. In fact, the 
fireball model for the emission of GRB predicts also the generation of 10 14 — 10 19 eU 
neutrino bursts. The Extreme Universe Space Observatory (EUSO), already ap- 
proved for accommodation study on the International Space Station would measure 
such UHEN, also in coincidence with the photons associated to the burst. As we can 
see, the near future might provide us with very good chances to detect such quan- 
tum gravity effects, or at least to set rather strict bounds on the theories predicting 
them. As a matter of fact, constraints upon the parameters defining such Lorentz 
violating theories have already been established by using current observations!! 
The extension of the modified dispersion relations (|l|) to other particles has pro- 
vided a way to circumvent some traditional astrophysical paradoxes such as the 
existence of the GZK cutoff, the observation of multi-TeV photons from Markar- 
ian 501 and the pion-stability paradox related to the structure of the air showers 
produced by high energy cosmic rays.tS The pioneering work revealing the appear- 
ance of Planck scale modifications to photon propagation in the framework of loop 
quantum gravity Hj was made by Gambini and PullinEj Also in the framework of 
loop quantum gravity the present author, in a collaboration with J. Alfaro and H. 
Morales-Tecotl, has extended such results by developing an heuristic semiclassical 
approach which allows the construction of effective field theories for photons and 
neutrinos, including Planck scale corrections arising as a manifestation of quantum 
gravitational effectsJiij In particular, dispersion relations of the form dl])_arise from 
such construction. An alternative approach inspired in string theory ta has been 
developed by Ellis et. al. and includes both photons and spin 1/2 particles.Ej The 
paper is organized as follows: in section 2 a brief review of the basic assumptions 
underlying the semiclassical approximation developed in Ref.(ll) is given. The re- 
maining sections contain new material which can be considered as the first steps 
to investigate a fully consistent coupling of the effective theories involved. Section 
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3 contains the formulation of the Gambini-Pullin electrodynamics with sources in 
terms of the standard electromagnetic potentials. In section 4, an action for the 
relativistic particle yielding Planck scale modified dispersion relations is presented. 
Finally, in section 5, 1 discuss the extension to a Dirac particle of the effective action 
previously found for two-components spin 1/2 fcrmions. 

2. The Effective Theories 

Each effective matter Hamiltonian is defined as the expectation value of the corre- 
sponding quantum gravity operator in a semiclassical mixed state which describes 
a flat metric together with the corresponding matter field. The requirements and 
properties of such a state are made precise in the sequel. The resulting effective 
theories violate Lorentz covariance at the Planck scale and such violation can be 
understood as a spontaneous symmetry breaking generated when taking the expec- 
tation value in the semiclassical state. 

In this section I summarize the procedure for the case of (the magnetic sector 
of) Maxwell theory. The starting point is the corresponding Hamiltonian 

H= [ d 3 x- -^=[E a E b + B a B\ (3) 
J s 2 Vdet<? 

where the space time is assumed to be a manifold M with topology £ x JR. Here 
£ is a Riemannian 3-manifold with metric q a b, El an d B_ a denote the electric 
and magnetic fields respectively. Thiemann has proposed a general regularization 
scheme that produces a sound mathematicaLdefinition for all the operators entering 
in the description of loop quantum gravity.ll3 Such regularized operators act upon 
states which are functions of generalized connections defined over graphs. A basis 
for such space is provided by the so called spin network states. A graph T is a 
set of vertices v £ V(T) in £ which are joined by edges e. The regularization 
procedure is based upon a triangulation of space which is adapted to each graph. 
This means that the space surrounding any vertex of T is filled with tetrahedra 
A having only one vertex in common with the graph (called the basepoint v(A)) 
plus segments s/(A) starting at A and directed along the edges of the graph. In 
the regions not including the vertices of T the choice of tetrahedra is arbitrary and 
the results are independent of it. The arcs connecting the end points of s/(A) and 
sj(A) are denoted by au(A) and the loop a/j := s/ o au o can be formed. 
A fundamental property of this procedure is the use of the volume operator V as 
a convenient regulator. In this way, the action of the operators is finite and gets 
concentrated only in the vertices of the graph. 

In the case of the magnetic sector of ([|) Thiemann's regularization leads to the 
operator 



e JKL e MNP x 
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x w lLA w lPA , (h ajK(A) - l) [h aMN{AI) - l) . (4) 

where 

w iL A = tr\Tih SL ( A ) 

and E(v) is a factor related to the number of edges that joint at the vertex v. In 
Eqs.([|[|) ^a(A) {ha(A)) denote parallel transport operators (holonomies) of either 
the gravitational connection Ai a or the electromagnetic connection A a , respectively, 
along the path a associated to the tetrahedron A. Such holonomies are SU(2) 
and U(l) group elements, respectively. The corresponding trajectories have been 
previously defined. 

Next, the quantum state that produces the semiclassical approximation is de- 
scribed. To this end let us consider an ensemble of graphs together with their 
adapted triangulation ( which means a set of segments {s/(A)} for each graph), 
characterized by some probability distribution P(T). To each graph T one asso- 
ciates a wave function \T,C,E_, B) which is peaked with respect to the classical 
electromagnetic field configuration together with a flat gravitational metric and a 
zero value for the gravitational connection. In other words, the contribution for 
each operator inside the expectation value can be estimated as 

(r,£,E,B\...q ab ...\T,£,E,B) 

{T,C,E,B\...A ia ...\T,C,E,B) 

while the expectation values including the electric and magnetic operators are esti- 
mated through their corresponding classical values E and B. Not surprisingly, the 
semiclassical state specifies both the classical coordinate and the classical momen- 
tum for each pair of canonical variables. The scale C » £p oi the wave function 
is such that the continuous flat metric approximation is appropriate for distances 
much larger that £, while the granular structure of spacetime becomes relevant 
when probing distances smaller that C. Such scale will have a natural realization 
according to each particular physical situation. In the sequel we set T = for 
simplicity and denote the semiclassical state as |T, S) — |T, C,E_,B). 

In a very schematic way I summarize now the method of calculation. For each 
graph r the effective Hamiltonian is defined as Hr = (r, S'l-ffrlT, S). For a given 
vextex, inside the expectation value, one expands each operator in powers of the 
segments s/(A) plus derivatives of the matter operators. In the case of (Q) this 
produces 

R r = E E (r,^|Z Pigi W...3 a ^.l M ( V )T ai ...^P^---( t ;, S (A))|r,5). 

vGV(r) v(A)=v 

(7) 




= Sab + O 

- »4 



c 

C 



(6) 
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where T contains gravitational operators together with contributions depending on 
the segments of the adapted triangulation in the particular graph. Next, space is 
considered to be divided into boxes, each centered at a given point x and with 
volume £ 3 as d 3 x. The choice of boxes is the same for all the graphs considered. 
Each box contains a large number of vertices of the semiclassical state (C » £p), 
but it is considered as infinitesimal in the scale where the space can be regarded 
as continuous. The sum over the vertices in (||) is subsequently split as the sum 
over the vertices in each box, plus the sum over boxes. Also, one assumes that 
the electromagnetic operators are slowly varying within a box (C « A, with A 
been the photon wavelength), in such a way that for all the vertices inside a given 
box one can write (r, S\ . . .F ab (u) . . . |T, S) = pF^^x). Here F ab is the classical 
electromagnetic field at the center of the box and p is a dimensionless constant 
which is determined in such a way that the standard classical result in the zeroth 
order approximation is recovered. Applying the procedure just described to (Q) 
leads to 

hf = J2 9i ( f ) • ■ ■ (° ai ■ ■ ■ -^-pM E £ p x 

Box w£Box 

]T p n+1 {T,S\^f ai ...^^-(v, S (A))\T,S). (8) 

v(A)=v P 

In the above, n + 1 is the total number of factors F pq (x) . The expectation value 
of the gravitational contribution is expected to be a rapidly varying function inside 
each box. Finally, the effective Hamiltonian is defined as an average over the graphs 
r, i.e. over adapted triangulations : H B = ^ r P(T)Hp. This effectively amounts 
to average the expectation values remaining in each box of the sum (||). We call 
this average T ai ,, pqpi qi — {x) and estimate it by demanding T to be constructed 
from the flat space tensors 5 a b and t a bc- In this way one is imposing isotropy and 
rotational invariance on our final result. Also the scalings given in ([|) together 
with the additional assumptions: (r, 5|...V'...|r, S) — ► £p,s° — ► Ip are used . 
Let us remark that the above average can be understood as taking the expectation 
value of the Hamiltonian (^) in a mixed state characterized by the density matrix 
p = J2r |r, S) P(r) (r, S\. After replacing the summation over boxes by the integral 
over space, the resulting Hamiltonian has the final form 

R B = J d 3 xJl piqi (x)... (V...^F p9 (x)) T a ,.™-(f). (9) 

Some comments are now in order. A rigorous semiclassical treatment of loop 
quantum gravity is still in the process of developmentO Since the approach pre- 
sented here has made use only of the main features that semiclassical states should 
have, all dimensionless coefficients in the expectation values that contribute to 
T ai ... pqpi 91 '"(f) in (||) remain undetermined. Besides, the calculation has not been 
performed in a covariant way. On the contrary, the results are expected to be valid 
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only in a preferred reference frame. Thus, the states considered so far will not 
annhilate the Hamiltonian constraint of quantum gravity. 

By applying the above procedure to photons (7) and to two-component massive 
neutrinos with definite chirality (u) the following results, including corrections to 



order £%, are obtained 113 



H 7 



d 6 x 



B 2 + E 2 ) + 6 8 £p ( B- (V x B) + E- (V x E) 



r,)i( 

+6> 3 Ip ( R a V 2 B a + E a V 2 £ a ) + 6 2 l P E a d a d b E b + 9 4 C 2 1\ ( jf ) 



W = / d A x 



l 7T(x) (l + Kl ^ + ...- K2 l p V 2 ^) T d d d i{x) + 

+ ^(£)^ + ...- K4 4vj^) 

+ jH T (x) (i* 2 ) (l + «5 Y + K 6 tp r a d a ) £(x) + c.c. 



(11) 

Here tt(x) = i£,*(x) and r a = —(i/2) a a with cr° being the standard Pauli matrices. 
In particular, the effective theories given in ( |Io| ) and (|ll]) imply Lorentz violating 
Planck scale modifications of the corresponding particle energy-momentum relations 
which are calculated in Ref.(ll). In both cases the scale C has been estimated as 
the De Broglie wavelength of the corresponding particle. In the case of neutrinos 
the condition m„ << p v was assumed. 

Before closing this section let us emphasize that the effective theories ([lC]) and 
( pd| ) are expected to be valid in a particular reference frame, the most natural one 
being that in which the CRB spectrum looks isotropic. This means that the in- 
volved scales lp and C will experience FitzGelrald-Lorentz contraction in going to 
the laboratory frame, for example. Also, the velocity of light will not have a uni- 
versal value, exhibiting corrections depending on ip which arise from the modified 
dispersion relations. An alternative point of view allowing for deformed dispersion 
relations which are^alid in every reference frame has been recently proposed^ and 
further elaborated] 2 !] This requires the formulation of a relativity principle having 
two observer independent scales: the speed of light constant and the Planck-length 
constant, which can be realized via non- linear realizations of the Lorentz group. An 
analysis of the common features and main differences between the approaches of 
Refs.(19) and (20) has also appeared.0 

3. The Gambini-Pullin Electrodynamics with Sources 

In this section I present the first steps towards a more detailed discussion of the 
modified electrodynamics obtained in (fol). In order to have the correct normaliza- 
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tion in the zeroth order case (£p — 0) it is convenient to make the field redefinition 

(1 + 6 7 (£ P /£) 2 ) 1/2 ^ — > Ei (12) 

and similarly for B^ — ► £?,;. Considering only the contribution linear in ip one is 
left with the effective Hamiltonian density 

U EM = \ (b 2 + E 2 ) + e s l P (B • (V x S) + E • (V x £)) , (13) 



2 

which was previously obtained by Gambini and PullinEl This theory predicts bire- 
fringence effects which are manifest through different propagation velocities for left 
and right polarized photons. Such velocity difference is proportional to the pa- 
rameter 8$ in Eq. ( [l3| ) and it is linear in £p. By analizing the presence of linear 
polarization in the optical and ultraviolet spectrum of some cosmological sources 
the bound 8 S < 10~ 3 is obtained.tJ 

Adding the appropriate sources in the first order action 

S[*, At, E 3 ] = J dtd 3 x (-Ei A, - H EM + $ {di Ei ~ 4tt P ) + AnJ.Ai) , (14) 

with B = V x A and the potential $ acting as a Lagrange multiplier, the modified 
Maxwell equations 

V • E = Anp, V • B, = 
Vx|b + 26> 8 £ P V x B) - ^ = An J, 

- / \ r) R 

Vx (E + 26 8 l P VxE) + — = 0, (15) 



are obtained. From Eqs. (|15|) one can prove the continuity equation for the electric 
charge as a signal of consistency. In order to write a second order Lagrangian 
formulation of Maxwell equations (JlJ) in terms of the basic fields A^ = (&,A) it is 
convenient to reintroduce the potentials starting, as usual, from the homogeneous 
equations. Here I use the definition F^ v = d^A v — d y A^ and the conventions of 
Jackson. The magnetic field retains the standard relation with the vector potential: 
Bi = CijkFjk, while the relation defining the electric field is changed to 

ldl 

E + 26 8 £ P Vx£ = -V$ — = F 0i . (16) 

c at 

Eq. ( filf ) implies that the resulting Lagrangian will be non local. In order to invert 
(Eq), the operator 



Mr.\ x , y ) = S^x - y) (% + 28 8 £ P e lkj d k ) = S(x° - y°)Mr.\x>, ff), (17) 
is defined in such a way that 

Ei{x) = f d 4 yM lJ (x 1 y)F Q] (y). (18) 
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The Lagrangian density is £ = — E l Ai — H, where the velocities Ai are introduced 
via the equation (|l8|), which defines as a non local function of the potentials. 
This leads to 

£ = - [e 2 - B 2 ) + 6$£ P (S ■ (V X E) - B ■ (V X B)j - 4tt J»A a . (19) 

The above Lagrangian density violates P, but preserves C and T. The CPT violation 
is produced because ( |l9| ) is both non local and non invariant under proper Lorentz 
transformations. Nevertheless, (|l^) is invariant under rotations and thus it is valid 
in a preferred coordinate system which we identify with that in which the CRB 
looks isotropic. It is interesting to observe that the Lagrangian density (|l|) is power 
counting non-renormalizable, thus constituting a case which is not considered in the 
general framework of Ref.(3) to discuss Lorentz and CPT violations. All effective 
theories generated via the framework presented here are expected to enjoy such non- 
renormalizability property. Also, ( |l9| ) describes an effective theory which should be 
valid only at energy scales much lower than the Planck mass. A convenient way of 
presenting the action arising from (119) is 



S = J <fxd% \ (Ei(x) M- q \x,y) E q (y) - B t (x) M~ 1 (x, y) B q {y))-^ J cfxJ^A, 

(20) 

where the operator M i(j 1 (x, y) satisfies the additional symmetry condition M iq 1 (x, y) 
= M~^(y,x). Varying the action ( pp[ ) leads indeed to the modified Maxwell equa- 
tions ([l5]). Before closing this section I discuss the explicit inversion of the operator 
Mjl (x, y). In momentum space, with the convention V = i k, one has 

(M~ x ) . m = S rm + (in£ P ) Eiim k t , c; = int Pl (21) 

which inverse is 

Mij = (Sij + seijpkp + <; 2 hkj) . (22) 

1 + q 2 k 2 

The pole in (|22| ) signals the need of a cutoff in order to regulate the Fourier trans- 
form. This is in accordance with the effective character of the theory, which is no 
more valid for momenta close to the pole position koo = 1/(k£p). In coordinate 
representation, (E2h leads to the Fourier transform 



M(z) = I d 3 k a** l — ^ = — fkdk - f n 2 i Z ,, 2 . (23) 

J 1 - n 2 l 2 P k 2 z J 1 - (k 2 £ 2 p ) k 2 

For energies which are low compared to fcoo, it will be enough to consider the local 
approximation to the operator Mjj (x, y) , which in turn leads to a simpler local 
effective theory. This is done by expanding the denominator of ( |2^ ) in power series 
of £p and integrating. The result, up to second order in £p is 

My (z) = (2tt) 3 S 3 (z) (<% + Kl P e ijp d p + k 2 £ 2 p (d.dj - S, tj V 2 ) + 0{k 3 £ p ) ) , (24) 
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with z = 



x-y. 



4. An Action for the Relativistic Particle 

In order to consistently couple the electrodynamics considered in section 3 to a point 
particle it is necessary to have an action for the latter which naturally incorporates 
corrections to the dispersion relations of the type discussed previously. A simple way 
to do this is by starting from a first order action where the corresponding modified 
energy-momentum relation is incorporated as a constraint. In other words, suppose 
that one needs to consider the modified dispersion relation 



F(p°)-(f+m 2 )=0. 



(25) 



The action to consider is 
S = 



A 



(26) 



where A is a Lagrange multiplier. The equations of motion are 
Sp : a 

8p l : -ii + Xp 1 = 



XdF 2*° 



P = A' 



(27) 



together with (25). Here G{u) — dF(u)/du. Substituting p° and p 1 back into the 
action yields 



dr ^G- 1 



2i° 
A 



2i° 
A 



A 



2A 



(28) 



The second order action is obtained by eliminating the auxiliary field A via its 
equation of motion arising from (^8|) . In order to illustrate the procedure in a more 
tractable situation I choose 



F(p°) 
G- 1 (y) 



(p°) 2 + a^(p°y 

1 



G(p° 



2p° 



a£pp° 



y, 



^(v/or^y-i 



In the approximation linear in £p I obtain 
G- 1 (y) 



12 



a£ P y 3 



(29) 



(30) 



Clearly, the dispersion relation arising from the choice (|29j) together with (|25|) re- 
produces Eq.(Q) in the zero mass limit. Substituting in ([28|), and after some algebra, 
one obtains 



dr 



v 

2A 



1 (x ) 3 



A 



ay 



(31) 
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The equation of motion for A is 



v 2 2 (±°y 



0, 



(32) 



which is also solved to first order in tp by making the ansatz A = Ao (1 +^pAi). 
The result is 

. 

= --a ■ , (33) 



An 



m 



1 m (i )" 

"3" (^2)3/2 



leading to the action 

S = mo I dt 



v/T_ ^~? (moM (i- i ^) 



(34) 



Observe that the limit lp = correctly reproduces the well known relativistic action 
for the point particle. From the previous equations one can write the energy and 
momentum in terms of the velocity as 



V 



™o7 - T7^ m o7 2 |l-;7 



m 7«'(l + ^KM7 3 ). (35) 



At this stage it is appropriate to compare the results (35) with those of Ref.(20), 
which, to first order in £p, are 



2„,2 



p = m 7 - £pm -f 



P 1 = rriQ^v 1 (1 - £pm j) . 



(36) 



The discrepancy is telling us that the action ( p4| ) is not an scalar under the specific 
non-linear representation of the Lorentz group proposed in Ref.(20). Finally, it is 
interesting to emphasize that Planck scale corrections to cither particle propaga- 



tion or quaii^L 
covariance.l 



Id interactions need not necessarily imply violations of Lorentz 



5. A Lagrangian for Dirac Particles 

The aim in this section is to construct a modified Lagrangian density for Dirac 
particles, starting from the theory given by (|TT|) . Such two component theory can 
be embedded into a four component realization by demanding '5 to be a Majorana 
spinor 









. x _ 




-i a 2 C _ 



In this notation, the equations arising from (|llj) are 
d 



(37) 



.... iAa-V 
at 

■ d . i „ 
i— +iA a • V 
at 



£-m(l-iClpa-V)x = 0, 
X-m(l-iC£ P a-V)£ = 0, 



(38) 



Quantum gravity corrections to particle ... 11 



with A = (l + Z3£pV 2 ) and C,D being constants. For simplicity I discuss here 
the rather unrealistic case C — ► oo. Next, x and £ are considered to be independent 
spinors. In the conventions where 75 is diagonal, (70) 2 = 1 and the signature is 
(H ), one can verify that 

(«7^+i(D4V 2 ) 7 • V -m(l-iCe P ^-v)) * = (39) 

reproduces the equations (j38|). The spin operator is given by E fe = (i/2)ekiml l l m ■ 



The Lagrangian density that yields Eq. (39) is 



C D = - (fytfJ + iD^ (V 2 *) 7 fc (9 fc *)-m§(l-iC^pS fe a fe )* + /i.c. 

(40) 

which is invariant under the global phase transformation (5^ = iSQ^. The asso- 
ciated Noether current, which will be identified with the electromagnetic current 
is 

J a = y 1 °y, J k = f( 7 fc + Cml P E fc )*, (41) 

to first order in lp. The coupling with the modified electrodynamics ( |l9| ) is made 
in the standard gauge invariant way via the replacement — > + ie^4 M in (|40|). 

Acknowledgments 

The author wishes to thank D.V. Ahluwalia and N. Dadhich for the warm hospi- 
tality extended to him in Pune. Also he acknowledges H. Morales-Tecotl for many 
illuminating discussions about loop quantum gravity. This work was supported in 
part by grants CONACyT 32431-E and DGAPA IN11700. 

References 

1. P. Huet and M. Peskin, Nucl. Phys. B434, 3 (1995); J. Ellis, J. Lopez, N.E. Mavromatos 
and D.V. Nanopoulos, Phys. Rev. D53, 3846 (1996). 

2. G. Amelino-Camelia, J. Ellis, N.E. Mavromatos, D.V. Nanopoulos and S. Sarkar, Nature 
393, 763 (1998). See also the contribution of S. Sarkar to these Procceedings. 

3. G. Amelino-Camelia, Nature 398, 216 (1999); Lect. Notes Phys. 541, 1 (2000). 

4. D.V. Ahluwalia, Nature 398, 199 (1999); GZ. Adunas, E. Rodriguez-Milla and D.V. 
Ahluwalia, Phys. Letts. B485, 215 (2000); ibid. Gen. Rel. Grav. 33, 183 (2001). 

5. D. Colladay and V.A. Kostelecky, Phys. Rev. D55, 6760 (1997), Phys. Rev. D58, 
116002 (1998) ;V.A. Kostelecky and CD. Lane, J. Math. Phys. 40, 6245 (1999);V.A. 
Kostelecky and R. Lehnert, Phys. Rev. D63, 065008 (2001); V.A. Kostelecky, arXiv. 



hep-ph/010422 r i and references therein; D. Colladay and P. McDonald, arXiv: hep- 
ph/020206$ ~ 



6. R. Bluhm, arXiv: hep-ph/0111325 and references therein. 

7. S.D. Biller et. al., Phys. Rev. Lett. 83, 2108 (1999). 

8. P.N. Bhat et. al., Nature 359, 217 (1992). 

9. V.A. Kostelecky and CD. Lane, Phys. Rev. D60, 116010 (1999); J.M. Carmona and 
J.L. Cortes, Phys. Letts. B494, 75 (2000); C. Lammerzhal and C. Borde, in Lecture 
Notes in Physics 562, Springer 2001 ; R. Brunstein, D. Eischler and S. Foffa, arXiv: 



12 Quantum gravity corrections to particle 



12 
13 



16 



17 



hep-ph/010630^ , S. Liberati, T. Jacobson and D Mattingly, arXiy : \hep-ph/01 10094 ; 
T. Jacobson, S. Liberati and D. Mattingly, arXiv: hep-ph/0 11220^ T.J. Konopka and 
S.A. Major, arXiv: hep-ph/ 020 1184; D. Sudarsky, L. Urrutia and H. Vucetich, arXiv: 



gr-qc/0204021. 



10. G. Amelino-Camelia and T. Piran, Phys. Rev. D64, 036005 (2001); G. Amelino- 
Camelia, P hys. Letts. B528, 181 (2002); J. Alfaro and G. Palma, arXiv: Iftep 



th/0111176 . 



11. For a review see for example C. Rovelli, in Living Reviews Vol 1, 1998-1, 
http://www.livingreviews.org/Articles. See also the contribution of G. Date to these 
Procceedings. 

R. Gambini and J. Pullin, Phys. Rev. D59, 124021 (1999). 

J. Alfaro, H. Morales-Tecotl and L. Urrutia, Phys. Rev. Lett. 84, 2318 (2000); J. Al- 
faro, H. Morales-Tecotl and L. Urrutia, Phys. Rev. D65, 103509 (2002); J. Alfaro, H. 
Morales-Tecotl and L. Urrutia in JHEP Proceedings, Cartagena de Indias 2000, High 
Energy Physics; J. Alfaro, H. Morales-Tecotl and L. Urrutia, in Proceedings of the Ninth 
Marcel Grossmann Meeting on General Relativity, eds. R.T Jantzen, V. Gurzadyan and 
R. Rufini (World Scientific, 2002). 

14. V.A. Kostelecky and S. Samuel, Phys. Rev. D39, 683 (1989); Phys. Rev. D40, 1886 
(1989); V.A. Kostelecky and R. Potting, Nucl. Phys. B 359,545 (1991), Phys. Lett. B 
381, 89 (1996) 

15. J. Ellis, N.E. Mavromatos and D.V. Nanopoulos, Gen. Rel. Grav. 32, 127 (2000); J. 
Ellis, N.E. Mavromatos and D.V. Nanopoulos, in Tegernsee 1999, Beyond the desert; J. 
Ellis, K. Farakos, N.E. Mavromatos, V.A. Mitsou and D.V. Nanopoulos, Astrophysical 
Jour. 535, 139 (2000); J. Ellis, N.E. Mavromatos, D.V. Nanopoulos and G. Volkov, 
Gen. Rel. Grav. 32, 1777 (2000). 

T. Thiemann, Class. Quan. Grav. 15, 1281 (1998); T. Thiemann, Class. Quant. Grav. 
15, 839 (1998); T. Thiemann, Phys. Letts. B380, 257 ( 1996). A comp lete and self- 
contained review of these ideas is in T. Thiemann, arXiv: gr-qc/01 10034- 
T. Thiemann, Class. Quant. Grav. 18, 2025 (2001); T. Thiemann and O. Win- 
kler, Class. Quant. Gr av. 18, 2561 (2 001); T. Thiemann and O. Winkler, arXiv: 



hey-th/0005 234 ; arXiv: hep-th/000523, 



qc/ 010203 f. 



H. Sahlmann and T. Thiemann, arXiv: 



gr- 



R.J. Gleiser and C.N. Kozameh, Phys. Rev. D64, 083007 (2001). 
G. Amelino-Camelia, Int. J. Mod. Phys. Dll, 35 ( 2002). 
J. Magueijo and L. Smolin, arXiv: hep-th/011209C. 



18. 
19. 
20. 

21. G. Amelino-Camelia, D. Benedetti and F. D'Andrea, arXiv: lfiep-th/0201245 . 

22. T. Padmanabhan, Phys. Rev. 57, 6206 (1998); K. Srinivasan, L. Sriramkumar and T. 
Padmanabhan, Phys. Rev. 58, 044009 (1998); S. Shankaranarayanan and T. Padman- 
abhan, Int. J. Mod. Phys. D10, 351 (2001). 



